Modeling of realistic cladding structures for air-core photonic band-gap fibers 
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Cladding structures of photonic band-gap fibers often have air-holes of non- 
circular shape and, typically, close-to-hexagonal air holes with curved corners 
are observed. We study photonic band-gaps in such structures by aid of a two- 
parameter representation of the size and curvature. For the fundamental band- 
gap we find that the band-gap edges (the intersections with the air line) shift 
toward shorter wavelengths when the air-filling fraction / is increased. The 
band-gap also broadens and the relative band-width increases exponentially 
with f 2 . Comparing to recent experiments [Nature 424, 657 (2003)] we find 
very good agreement. 
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In air-silica photonic crystal fibers (PCFs) an arrange- 
ment of air-holes running along the full length of the 
fiber provides the confinement and guidance of light in 
a defect region. For photonic band-gap (PBG) guiding 
PCFs the air-holes have been arranged in various ways 
such as in a triangular latticei, but honey-comb^ and 
kagome 3 arrangements are other options. Cregan et al^ 
have demonstrated that light can be guided in an air- 
core by means of the PBG effect and this observation 
has stimulated an avalanche of both basic and applied 
research. For recent reviews we refer to Refs. I ll5l and 
references therein. 

From a modeling point of view the air-holes are gen- 
erally assumed circular^i and for fibers fabricated by 
the stack-and-pull method this is also typically what 
is observed for moderate air-filling fractions. 8 However, 
for the high air-filling fractions / > 0.8 recently em- 
ployed for PBG guidance the air-holes tend to be non- 
circular Aifl From a simple geometric consideration it is 




Fig. 1. The right panel shows a scanning-electron mi- 
crograph of the cladding structure of a PBG fiber with 
a pitch of A ~ 2.7 ^m. The left panel illustrates the 
two-parameter representation of the air-hole shape. 



easily shown that the theoretical upper limit of the air- 
filling fraction in a structure of close packed circular voids 
is / = 7r/(2v3) — 0.91. Hence, design and fabrication 
of air-core PBG fibers call for a detailed modeling of the 
spectral position of the PBGs which goes beyond the as- 
sumption of circular air holes. In this Letter we consider 
the triangular cladding arrangement (with a pitch A) first 
employed by Cregan et ai.ji but for the air-holes we take 
the non-circular shape observed recently in Refs. Ij^lTol into 
account. From scanning-electron micrographs (SEMs) 
we find that the air-holes are of over-all hexagonal shape 
with curved corners, see Fig. ^ (this conclusion was also 
emphasized in Ref . llUS) . This suggests that the holes can 
be parametrized by two parameters; the edge-to-edge dis- 
tance D (corresponding to the diameter of a circular air- 
hole) and the diameter of curvature, D c , at the corners 
(see left panel of Fig. ^| . In this two-parameter represen- 
tation the air-filling fraction is given by 
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where Ah ex (x) = v / 3a; 2 /2 and A c { TC {x) — ir(x/2) 2 are the 
areas of a hexagon (with edge-to-edge distance x) and a 
circle (of diameter x), respectively. 

For the optical properties we apply Maxwell's equa- 
tions to a macroscopic and isotropic loss-less dielectric 
medium and assume a linear dependence of the displace- 
ment field on the electrical field. We consider a harmonic 
mode H(r, t) — J? w (r)e iwt with angular frequency ui and 
substituting into Maxwell's equations the magnetic-field 
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Fig. 2. Plot of (3/k versus k± along principal direc- 
tions in the Brillouin zone for (3A = 14.45 and a struc- 
ture with D/A = 0.96 and D c /A = 0.55 (see inset of 
Fig-EJ- The filled curve shows the corresponding density- 
of-states (the projection of the data onto the y-axis). 

vector is then governed by the wave equation^ 

Vxi-Vxff u (r) = ^H u (r). (2) 
e(r) c z 

Here, c is the velocity of light and e is the dielectric func- 
tion which we in the following assume independent of 
frequency by which the wave equation becomes scale- 
invariant. All the results to be presented can thus be 
scaled to the desired value of A. 

For a PCF geometry, the cladding dielectric function 
e(r) is periodic in the transverse plane and translational 
invariant along the fiber-axis (the z-axis) . The solution is 
according to Bloch's theorem a plane wave modulated by 
a function h^lr^) with the periodicity of the dielectric 
structure in the transverse direction 

H u (r) = h u) (r±)exp(ik±-rx + i^z). (3) 

Substituting this ansatz into Eq. © we get an eigenvalue 
problem for u>(kj_ , (3) which we solve by the aid of a plane- 
wave basis (typically 128 x 128 plane waves) with periodic 
boundary conditions. 12 For the dielectric function we use 
e = 1 in air and s — (1.444) 2 = 2.085 in silica. From a 
computational point of view our simulations thus follow 
the same lines as those used for structures with circular 
holes£ 

In Fig. [2] we show photonic bands calculated for (3 A = 
14.45 and a structure with D/A = 0.96 and D c /A = 0.55 
(see inset of Fig.[3Jl. It is common to introduce the free- 
space wave- number k = 2ir/\ = uj/c and plot the ef- 
fective index (3/k — c(3/uj rather than the angular fre- 
quency uj. The bands are plotted as a function of k± 
along the usual principal directions in the Brillouin zone; 
from the T-point to the M-point and via the if-point 



back to the T-pointpi As seen the structure exhibits a 
band-gap which is particular clear from the filled curve 
which shows the corresponding density-of-states (the pro- 
jection of the data onto the y-axis). In this example 
the band-gap is around the index of air {(3/k — 1), but 
in general the position and width of the band-gaps de- 
pend on the propagation constant (3. In Fig. [31 we show 
the results of changing (3. The data-points result from 
a sampling of the Brillouin zone along the above men- 
tioned principal directions. The PBGs are those regions 
with no data-points and the solid lines indicate the cor- 
responding band-gap edges. The dashed line shows the 
air-line which passes through several band-gaps. Usu- 
ally the band-gap at the longest wavelength (the lowest 
frequency) is refereed to as the fundamental band-gap 
and the other gaps are denoted higher-order band-gapsr 6 
The slopes of the band-gaps are relatively large which 
suggest that e.g. scattering and bending-loss will in- 
crease almost abruptly when the wavelength approaches 
the band-edges. 

For PBG fibers with an air-coreii2ii£ the band-edges 
of guided modes will to a good approximation be given 
by the intersection of the air-line with the cladding band- 
edges in Fig. |3 For the fundamental band-gap we denote 
the upper and lower intersections by X u and A/, respec- 
tively. We have calculated the band edges for various 
values of D and D c ; Fig. ^summarizes the results. The 
band-diagrams of course depend on the two parameters 
D and D c , but we find that a single-parameter descrip- 
tion of the band-edges is possible in terms of the air-filling 




Fig. 3. Plot of (3/k versus A/A for a structure with 
D/A = 0.96 and D c /A = 0.55 (see inset). The data- 
points result from a sampling of the Brillouin zone along 
the usual path from the T-point to the M-point and via 
the if-point back to the T-point. The solid lines indicate 
the band-gap boundaries and the dashed line (3/k = 1 
shows the air line which passes through several band- 
gaps. 
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Fig. 4. Plot of the band-gap edges X u and A; of the fun- 
damental band-gap as a function of air-filling fraction /. 
The dashed lines are guides to the eyes. Recent experi- 
mental data™ is also included. 



Fig. 5. Relative band- width versus / 2 . The straight 
dashed line indicates the close-to-exponential depen- 
dence, Eq. 0. Recent experimental data— i is also in- 
cluded. 



fraction /, Eq. 

When the air-filling fraction is increased the center 
A c = (A u + A;)/2 of the band-gap shifts toward shorter 
wavelengths. Furthermore, the band-width AA = A u — A; 
at the same time increases and the relative band-width 
shows a close to exponential increase (see Fig.EJ, i.e. 

AA/A C ~ £ x cxp(.4 x / 2 ), (4) 

where A and B are positive numerical coefficients. The 
results in Figs. 14151 support the choice of high air- filling 
fractions for practical fibers and as seen we find very good 
agreement between our numerics and recent experimental 
results by Smith et al*& 

In summary, we have shown that realistic PBG 
cladding structures^iiS can be represented by an "ideal" 
two-parameter description which facilitate detailed nu- 
merical modeling. For the fundamental band-gap the 
band-gap edges (the intersections with the air line) shift 
toward shorter wavelengths for an increasing air-filling 
fraction / and the band-gap also broadens significantly. 
This observation may make air-core PBG fibers realistic 
for wavelengths even shorter than the 850 nm reported 
recently by Bouwmans et aV 
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